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A multispeckle technique for efficiently measuring correctly ensemble-averaged intensity autocor- 
relation functions of scattered light from non-ergodic and/or non-stationary systems is described. 
The method employs a CCD camera as a multispeckle light detector and a computer-based cor- 
relator, and permits the simultaneous calculation of up to 500 correlation functions, where each 
correlation function is started at a different time. The correlation functions are calculated in real 
time and are referenced to a unique starting time. The multispeckle nature of the CCD camera 
detector means that a true ensemble average is calculated; no time averaging is necessary. The tech- 
nique thus provides a "snapshot" of the dynamics, making it particularly useful for non-stationary 
systems where the dynamics are changing with time. Delay times spanning the range from 1 ms to 
1000 s are readily achieved with this method. The technique is demonstrated in the multiple scat- 
tering limit where diffusing-wave spectroscopy theory applies. The technique can also be combined 
with a recently-developed two-cell technique that can measure faster decay times. The combined 
technique can measure delay times from 10 ns to 1000 s. The method is peculiarly well suited for 
studying aging processes in soft glassy materials, which exhibit both short and long relaxation times, 
non-ergodic dynamics, and slowly-evolving transient behavior. 

PACS numbers: 07.60.-j, 78.35.+C, 81.40.Cd 



I. INTRODUCTION 



Dynamic light scattering has proven to be a powerful 
tool for investigating the dynamics of colloidal suspen- 
sions, emulsions, foams, and other complex fluids. Over 
the past decade, dynamic light scattering (DLS) has been 
increasingly applied to glassy systems with slow, and of- 
ten non-ergodic, dynamics. Slow dynamics pose experi- 
mental difficulties for DLS, particularly with respect to 
the stability of the lasers, and sometimes for the light 
detection system. Nonergodic samples pose even more 
serious experimental difficulties, which are exacerbated 
by the presence of slow dynamics. One way of address- 
ing the problem of long data acquisition times is to use 
an array detector. More importantly, an array detector 
can also address the problem of non-ergodic samples if 
the optics are designed so that the different elements in 
the array detector constitute essentially statistically in- 
dependent samplings of the system under study. 

In this article, we describe a method which employs 
a charge-coupled device (CCD) array detector for per- 
forming dynamic light scattering on systems that ex- 
hibit strong multiple light scattering. While the basic 
method we describe is applicable to single as well as 
multiple light scattering, our focus here is on multiple 
dynamic light scattering, also known as diffusing-wave 
spectroscopy (DWS) (I], 0, 0], where the methods we de- 
scribe are most effective. 



II. BACKGROUND 

In conventional DWS (or DLS), the basic idea is to 
illuminate a sample with a coherent beam of light from 
a laser and to measure the temporal fluctuations in the 
resulting speckle pattern of scattered light. In a typical 
experiment, the optics for collecting the scattered light 
is arranged in such a way that the intensity of a single 
speckle is measured. By focusing on a single speckle, the 
amplitude of the fluctuating intensity is optimized to ob- 
tain a strongly fluctuating signal. The temporal duration 
of a typical fluctuation, i.e. the lifetime of a speckle, is 
determined by the time it takes a particle to move such 
that the phase of the scattered light is changed by ap- 
proximately unity (or it). The dynamics of the scattered 
light, and hence the dynamics of the scatterers, can be 
more quantitatively characterized by the temporal auto- 
correlation function of the scattered intensity. In many 
cases, the autocorrelation function can be analyzed so 
that the time dependent mean square displacement of 
the scatterers is extracted ||. 

There are two important keys to the success of the con- 
ventional analysis that we wish to focus on here. First, 
the process must be stationary, at least for the duration 
of the experiment. That is, the dynamics of the system 
cannot change significantly during this time period. For 
example, the gelation process of a colloidal suspension 
can be monitored provided that the particle dynamics 
slow down (due to the increase of the elastic modulus) oc- 
curs on time scales longer than the measurement period. 
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In general an accurate measurement of the correlation 
function requires that the fluctuations in the intensity 
be measured over many fluctuation periods (e.g. ~ I0 4 
speckle lifetimes to achieve ~ 1% accuracy). Thus, in this 
case, the particle dynamics must not change significantly 
over a time period of ~ 10 4 speckle lifetimes. 

Second, the system must be ergodic. For a light scat- 
tering experiment, a system is ergodic if each speckle 
samples the same distribution of light intensities as the 
entire ensemble of speckles over the period of time that 
data is acquired (i.e. over the duration of the experi- 
ment). While this condition is satisfied for a large num- 
ber of systems, there are certain systems with very slow 
glassy dynamics that are of interest but which, neverthe- 
less, are not ergodic in the sense described above. Such 
systems pose a problem for light scattering and special 
steps must be taken in order to assure that meaningful 
correlation functions, which represent dynamics charac- 
teristic of the entire system, are measured ^, ||, ||] . 

The problem for light scattering presented by noner- 
godic samples was first understood in a clear fashion by 
Pusey and van Megen who, at the same time, presented 
a clever method of obtaining correctly averaged temporal 
autocorrelation functions 0] . This technique has proven 
to be very useful but has some drawbacks: it sacrifices 
some signal-to-noise in order to obtain a useful distribu- 
tion of light intensities, and the dynamics of the system 
must not change significantly over the duration of the 
experiment [Q, [^, [|. Moreover, it suffers from the same 
limitation as normal DLS experiments in that data must 
be collected over a time scale that is much longer that 
the longest time scale of the fluctuations one is trying to 
measure. 

Other schemes have been developed which overcome 
some of these difficulties, but not without a significant 
cost. The sample can be translated or rotated, for ex- 
ample, so that an ensemble average is obtained by brute 
force M. However, the rotation or translation introduces 
an artificial characteristic time scale into the correlation 
function, namely the time scale over which the system 
samples new speckles by virtue of the motion. Care must 
be taken to assure that this time scale is much longer than 
any dynamical time scale of the system one is interested 
in measuring. And once again, the duration of the ex- 
periment must be such that averaging occurs over a time 
scale that is long compared to the time scale introduced 
by the sample motion (this amounts to assuring that a 
sufficiently large number of speckles are sampled - e.g. 
~ 10 4 - to obtain the desired level of accuracy - e.g. 
~ 1%). Another method, developed recently by Schef- 
fold et at. @, H ||], uses two cells, one containing the 
sample with the slow nonergodic dynamics that is being 
studied, and another containing diffusing particles with 
ergodic dynamics, which serves to spatially average the 
nonergodic signal. This technique is only useful for ex- 
periments involving strong multiple scattering, but nev- 
ertheless provides a simple and efficient means of dealing 
with nonergodic systems. However, the only time scales 



probed in the primary cell are those that are shorter than 
the diffusive time scale set by the secondary cell. 

An attractive alternative to the techniques described 
above for dealing with nonergodic samples is to use a 
CCD array, conveniently contained in a commercial CCD 
camera, as a light detector 0, [Q[ Q [jj, 0. This 
allows the intensity from a large number of speckles to 
be simultaneously and independently measured. A time 
series of images of scattered light obtained from the CCD 
camera are stored in a computer and processed such that 
each pixel in the CCD array serves as an input to its 
own correlator. Since the fluctuation of each speckle is 
independent of the other speckles, processing scattered 
intensity data collected with a CCD camera is like have 
N independent samplings of the data, where N is the 
number of speckles in the field of the CCD array. If N 
is sufficiently large, e.g. N ~ I0 4 , then a complete en- 
semble average is obtained without any need to sample 
the system over times any longer than the actual dy- 
namical time scales that are under study. Thus, a direct 
ensemble average can be obtained without any need for 
time-averaging whatsoever. This approach solves several 
problems at once. First, it solves the nonergodicity prob- 
lem because a full ensemble average (of ~ I0 4 samples) 
is obtained. Second, the process under study need not be 
a stationary process because no time averaging is neces- 
sary, provided a sufficiently large number of speckles are 
used. Thus, an essentially instantaneous snapshot of the 
system dynamics can be taken. If those dynamics change 
in time, they can be followed simply by measuring tempo- 
ral autocorrelation functions at different starting times. 

There are two practical limitations of the CCD- 
camera-based technique. First, the minimum delay time 
for which the correlation function can be measured is set 
by the frame rate of the camera. In the experiments de- 
scribed in this article, we use a frame rate of 500 frames 
per second. Thus, delay times as small as a few millisec- 
onds are accessible. It is interesting and important to 
note that real time processing of the images is necessary 
since, to our knowledge, no commercial buffer memory is 
able to store them for the entire sequence of frames ob- 
tained in the course of a single experiment, which is typ- 
ically on the order of 10 5 to I0 7 frames. Recent advances 
in microprocessor speed allows real time computation of 
the correlation function. Hence the limiting factor for 
frame rate is the PCI bus velocity for transferring data 
from the acquisition board to the computer RAM. 

The second limitation of the technique is the limited 
sensitivity of the CCD array used, as compared to that 
of a photomultiplicr. This means that relatively strong 
scattered light signals are needed. In practice a power 
of 400 mW gives a satisfactory signal to noise ratio for a 
frame rate of 500 per seconds. 

In this article we first describe the CCD optical set up 
and calculational scheme for our computer-based corre- 
lator. In the next section we explain how a single camera 
can be used to take up to 500 correlation functions si- 
multaneously, with the reference starting time shifted by 



different amounts. We then show how the CCD camera 
technique can be coupled with the double cell technique 
to extend the time window over which measurements can 
be made from 10~ 8 s to 10 4 s. 



III. EXPERIMENTAL SETUP 

The sample is illuminated using polarized light from an 
argon-ion laser operating at a wavelength of A = 514 nm. 
Its maximal output power is 1.6 Watts. The laser beam 
is expanded to a diameter of approximately 1 cm and is 
incident on the sample cell, as shown in Fig. [l]. Multi- 
ply scattered light is collected by a 50-mm Nikon cam- 
era lens and focused onto a iris diaphragm. The lens is 
set up so that a one-to-one image of the scattered light 
from the output plane of the sample cell is formed on 
the diaphragm. The imaging of the output plane of the 
sample onto the pinhole provides a convenient means for 
collecting data from different parts of the sample. This 
feature is particularly useful in samples like gels or very 
viscous solutions where defects or small particles of dust 
may hinder the measurement; this set-up allows one to 
avoid defected regions of the sample. A polarizing cube 
is placed between the lens and the diaphragm to protect 
the camera from stray light. 

The low sensitivity of CCD arrays relative to photo- 
multiplier tubes means that the sample must be strongly 
illuminated. However, the incident light intensity must 
not be so high as to heat the sample. Therefore, it is 
crucial to collect as much scattered light as possible. We 
use a camera lens of 50 mm with a numerical aperture of 
1.7. 

The scattered light is detected by a Dalsa CCD cam- 
era, model CAD1-128A, which is placed approximately 
d ~ 15 cm behind the diaphragm. It is an 8-bit cam- 
era which we run at 500 frames per second. The images 
are transferred to a computer running at 500 MHz us- 
ing a National Instruments data acquisition board, model 
PCI1422. The combined speeds of the computer and data 
acquisition boards are sufficiently high to allow data to 
be analyzed in real time. 

Data analysis includes calculation of the ensemble- 
averaged temporal correlation function of the scattering 
intensity J, 
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and the ensemble-averaged spatial correlation function 

(J(r) J(r + dr)> 



C(Sr) 
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(2) 



where to is the time when the reference intensity is mea- 
sured, t > t ) is the delay time, r is the coordinate of 
a pixel on the detector, and r + <5r is the coordinate of 
another pixel a distance <5r away.j^, represent the average 
over the CCD chip. The denominator of the expression 
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FIG. 1: Experimental set up for the transmission geometry. 
The laser wavelength is A = 514 nm. We use a Dalsa Camera 
running at 500 frames/sec and a 500-MHz PC. 



on the right hand side of Eq. (|l|) is the square of the av- 
erage intensity determined according to the method sug- 
gested by Schatzel |ll| , in order to obtain a slightly more 
accurate measurement of <?2(Mo)- The algorithms used 
to calculate these correlation functions are described in 
greater detail below. 



A. Speckle size adjustment 

The size s of the speckles of scattered light formed on 
the CCD chip is given by 



s ~ dX/a 



(3) 



where d is the distance between the diaphragm and the 
CCD chip, A is the wavelength of light, and a is the 
diameter of the diaphragm. By adjusting d and/or a, the 
size of the speckles relative to the size of a pixel on the 
CCD chip can be controlled. We use a 2 mm x 2 mm 
CCD chip with 128 x 128 pixels. Each pixel occupies 
approximately 256 /im 2 with a 100% fill factor. 

To maximize the signal-to-noise, we need to collect 
data from as many speckles of scattered light as pos- 
sible. However, if the speckles are too small, for example 
if several speckles occupy a single pixel, optical contrast 
is lost and the signal for determining a temporal auto- 
correlation function gi (t) is degraded. Therefore, it is 
important to adjust the speckle size so that we obtain 
a satisfactory optical contrast with as much speckles as 
possible analyzed. 

In order to determine the number of speckles per pixel, 
we adjust the collection optics (i.e. the size of the iris) and 
measure the spatial intensity autocorrelation function 



C(p) 



(4) 



where Ii is the intensity of scattered light at pixel i and 
Ii+p is the intensity of light p pixels away. The corre- 
lation function C(p) is determined by averaging over all 
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FIG. 2: Spatial correlation function C(p) obtained using dif- 
ferent speckle sizes. The first minimum of this function is 
taken as the speckle spot diameter. 

pixels. Figure || shows C(p) for three different sets of 
light collections optics. Note that C(p) first decreases 
and then oscillates about a baseline intensity. We arbi- 
trarily define the position of the first minimum to be the 
average speckle size. Hence the speckle size is defined 
as the distance over which the intensity profile remains 
spatially correlated. 

To determine optimal configuration of the collection 
optics, we use the static scattering pattern from an im- 
mobile glass frit. The frit is sufficiently opaque that 
diffusion approximation for light transport through the 
medium is valid. The intercept for the intensity corre- 
lation function (i.e. 32(0)) is expected to lie between 1 
and 2. In order to quantitatively evaluate its value, a 
series of 2000 intercepts was measured. Measurements 
were performed by self-correlating 2000 images from the 
glass frit. For each image the mean intensity was calcu- 
lated as well as the mean squared intensity. The ratio of 
theses two quantities corresponds to 32(0) for this image. 
Six series of such measurements were made, each for a 
different speckle size. Figure ||(a) shows the six distribu- 
tions of 32(0). The average value (32(0)) for each speckle 
size is an evaluation of the contrast, whereas the spread 
A<7 2 (0) indicates how quickly (i.e. how many images must 
be correlated before) 32 (0) converges to its proper statis- 
tical value. As such, A32(0) quantifies the statistical and 
instrumental noise. To determine the best compromise 
between obtaining the highest signal with the lease noise, 
we plot (32(0) — 1)/A<7(0) as a function of speckle size. 
Figure ||(b) shows that the best configuration is obtained 
for a speckle size of 3 pixels in diameter. All subsequent 
measurements are performed in this configuration. The 
average intercept (32(0)} allows the determination of the 
normalization factor given for the Siegert relation: 

3 2 (i) -1 = /%1 (*)| 2 (5) 
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FIG. 3: (a) Histogram of the intercepts (<72(t) = 0) for 
2000 correlation functions for a frit at different speckle size. 
The theoretical intercept should be 2. (b) "Signal-to-noise" 
(172(0) — 1}/A(?(0) as a function of speckle size where (32(0)} 
is the average value of 32(0) and A$(0) its spread. It shows 
that the best compromise between lost of statistics and lost 
of contrast is obtained for a speckle diameter of 3 pixels. 

malized electric field autocorrelation function. We take 
P = 0.82 ± 0.05 in subsequent measurements and display 
our results in terms of the squared electric field correla- 
tion function |3i(i)| 2 using the Siegert relation above. 

IV. CALCULATION ALGORITHM 

In this section we explain how the bitmapped image is 
processed to calculate the autocorrelation function. 

A. Correlation Function Calculation 

At time t, a single CCD frame of scattered light is saved 
in the computer memory (RAM) . We call /™™ (t) the raw 
intensity of the light incident on pixel i for the image 
taken at time t, and subtract off an average background 
(ij- ) for pixel i due to stray light and dark counts: 



where 31(f) = (E(t )E(t + t)) / (\E(t )\ 2 ) is the nor- 



h(t) = r i aw {t) - (i* ack ) 



(6) 
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FIG. 4: Comparison of |<?i(t)| 2 obtained from a sample of 
0.5 /xm-diameter beads in pure glycerol for standard single- 
correlator DWS set up (— ) and multi-speckle camera based 
set up (Square). The inert shows the correlation function 
with the baseline subtracted demonstrating that reliable data 
is obtained down to 10~ 2 . 



We measure (lj? ack } by averaging 1000 images for each 
pixel with the laser turned off. Further refinements tak- 
ing into account other sources of extraneous signal can 
be found in ||l3f . While such refinements prove useful for 
DLS in the single scattering limit, we found that such 
corrections to be insignificant for DWS. 

One practical problem encountered in acquiring data 
is that some images may be skipped by the acquisition 
card due to synchronization problems. Thus, as a pre- 
caution, we always record the time t at which each image 
is obtained. 

Applying Eq. @ to our data, we calculate the corre- 
lation function: 



92 



'(*o,*) 



(Ji(* )/i(* + *o))i 



(7) 



Note that the average is performed over pixels rather 
than using the customary method of averaging over 
the initial time to- That is, the correlation function 

92 : (*Qj*) is calculated by multiplying two frames (one 
at time to , the other at to + *) together pixel by pixel, 
summing the results from all pixels, and then dividing by 
the number of pixels. In general, one frame suffices to 
obtain a good average for <72™^ (*<)>*) because it consists 
of a large number of pixels, 16384 = 128 2 in this case, 
which constitutes several thousand independent samples 
from which to determine g\ (to, t). 

The correlation function g^ (*o , *) is obtained in this 
manner is obtained from an ensemble average of the 
sample, with no need for the time averaging method 



typically used with conventional light detection. When 
data is taken from an ergodic system, both methods give 
the same result, as illustrated in Fig. ^ where we com- 
pare correlation functions obtained from a CCD camera 
to one taken with a conventional single detector-single 
correlator for a sample of polystyrene beads in glycerol 
obtained undergoing Brownian motion. Significant ad- 
vantages accrue for the CCD technique, however, when 
measurements are performed on a non-ergodic sample. 
Moreover, the CCD method can provide important ad- 
vantages when measuring correlation functions for non- 
ergodic systems. 



B. Multiple correlation function calculation 

In the preceding section, we discussed how a two-time 
correlation function 52 (to, *) can be obtained, to being the 
time at which the reference image is acquired and to + 1 
being a later time at which another image is collected and 
correlated with the reference. Such a correlation function 
provides a snapshot of the dynamical state of the system 
at a particular starting time, namely to, and is possible 
only because sufficient ensemble averaging is achieved by 
correlating thousands of pixels between two frames from 
a CCD camera. The ability to obtain such snapshots is 
particularly useful for studying systems whose dynamics 
are time-dependent, i.e. non-stationary systems. Exam- 
ples include systems which are gelling, coarsening, aging, 
or undergoing some other non-equilibrium non-stationary 
process. By measuring 52 (to, t) as a function of the start- 
ing time to, the evolving dynamics of such systems can 
be followed and studied. 

Here, we explain a scheme for following the evolving 
dynamics of a system by measuring many different corre- 
lation functions in parallel, each with a different starting 
time to- A unique and important feature of this scheme is 
that new correlation functions can be started at intervals 
as short as 1/500 s, while still collecting frames for cor- 
relation functions started at earlier times. In fact, up to 
500 correlation functions can be collected simultaneously, 
each with a different starting time to, and without any re- 
quirement that a previously started correlation function 
be finished before starting a new one. 

The key to the success of the scheme is the realization 
that in most cases, it is not necessary to measure the cor- 
relation function with the same temporal resolution for 
all delay times. Typically, small time steps are required 
for short delay times while larger time steps are required 
for longer delay times. The scheme we employ takes ad- 
vantage of this fact. It is essentially an interlacing scheme 
in which successive frames are sometimes used in differ- 
ent correlation functions, as will be explained below. In 
effect, the scheme represents a compromise between mak- 
ing maximal use of the information available in the mea- 
surements for calculating the correlation functions and 
obtaining and processing the data from the CCD camera 
frames at a rapid rate. 
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This compromise is necessary because calculating the 
correlation function between two frames is time consum- 
ing; only one pair of frames can be processed (multiplied 
and averaged) between the acquisition of two successive 
images. Image stacking is not practical because it quickly 
leads to memory overflow. Therefore, we employ an algo- 
rithm for calculating many correlation functions, labeled 
by j = 1, 2, 3, ... , jmaxi simultaneously in real time, that 
makes only the following manageable demands on the 
computer: 

• only one frame multiplication is performed between 
the acquisition of successive frames, 

• storage of a maximum of j max < 500 reference im- 
age frames plus the most recently acquired frame 
in memory, and 

• storage of a maximum of jmax correlation functions 
in memory, each with on the order of 100 channels. 

For a frame rate of 500 frames/s from a 128 x 128 pixel 
CCD camera, these demands can be met by a computer 
with a Pentium processor running at 500 MHz and pos- 
sessing 512 MB of RAM. An important feature of the al- 
gorithm is that it requires that only the reference frame 
for each correlation function, along with most recently 



Note that the spacing between the points in a correlation 
function is constant over the first N points but increases 
quadratically at later times. The time interval between 
two successive correlation functions j and j + 1 is ap- 
proximately N / f where / is the frame rate. This time 
interval can be tuned in order to best accommodate the 
rate of evolution of the system under study. For N = 1 
this time interval can be as short as the inverse frame 
rate. 

In Fig. |6| we show a series of correlation functions ob- 
tained using the multispeckle CCD camera technique de- 
scribed above. The correlation functions were obtained 
from a sample consisting of 0.5 wt.% alumina (AI2O3) 
beads, approximately 160 nm in diameter (AKP30 Sum- 
itomo), that were embedded in a crosslinked elastic net- 
work of polybutylacrylate. The first correlation function 
was started at a time to — immediately after a step 
strain of 5% was applied to the sample. Thus, the corre- 
lation functions measure the subsequent random motion 
of the beads that arises as a response to the stress within 
the sample. As the stress relaxes over the course of tens 
of minutes, the random motion of the beads becomes 
slower. This slowing down is reflected in the increased 



acquired frame, be stored in the computer at any mo- 
ment in time; this keeps the number of frames that must 
be stored in memory to a manageable number. This also 
allows all correlation functions to be calculated in real 
time as the experiment proceeds. 

The scheme we use is illustrated graphically in Fig. 
||. The vertical axis represents the different correlation 
functions labeled by j = 1, 2, 3, . . . , j ma x- The horizontal 
axis represents the successive frames that are acquired, 
labeled n = 1, 2,3, . . . with each frame being nominally 
separated in time by the reciprocal 1// of the frame rate. 
The reference image for the j th correlation function is 
given by image number 

n=n[ j) = l + (A + l)(j-l) . (8) 

where the subscript "1" indicates that this is the first 
frame of the j th correlation function. Each of these ref- 
erence images is stored in memory; N is typically of order 
10. The reference image and the next N — 1 images in 
sequence are correlated with the reference image. The 
remaining frames used to calculate the j th consists of a 
sequence of selected images as illustrated in Fig. ||. The 
sequence numbers of the frames used to calculate the j th 
correlation function is summarized by the formula: 



(9) 

I 

decorrelation time of correlation functions started at in- 
creasingly later times after the initial step strain. In Fig. 
we plot the decorrelation time T1/2, defined as the time 
it takes a correlation function to decay to half its initial 
value, vs. time, referenced to the starting time to- It ex- 
hibits a power law behavior over two decades. A detailed 
discussion of this experiment will be published elsewhere. 



V. FAST AND SLOW TIME SCALES 

The shortest decorrelation time T1/2 measured in the 
experiments described in the previous section is fairly 
slow, approximately 0.1s. In this particular instance, the 
frame rate of the camera is sufficient to capture all the 
dynamics of the process in the correlation function. This 
is not always the case, however. Many materials, includ- 
ing polymer networks, colloidal gels, and glassy materi- 
als possess two characteristic ranges of relaxation times, 
where only the slower of the two is sufficiently slow for 
the multispeckle CCD camera method to work. In these 
cases, the faster relaxation processes are not measured. 

To illustrate such a process, we use the multispeckle 



f n[ j) +i-l = (N+l)(j 
1 (j - 1) + - 1) - ±N{N - 1) 



if 1 < i < N 
if i > N 
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FIG. 5: Schematic representation of image selection and use in the multicorrelator scheme. For this illustration N — 5. The 
horizontal axis is the image number n (or n' J ') and the vertical axis is the correlator number j. The numbers inside the boxes 
indicate the i th image used in correlator j. 
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FIG. 6: Correlation functions |gri(io,£)| 2 taken after a shear 
step of 5%. Multiple correlation functions are obtained in real 
time with the ones started later having larger decorrelation 
times; the decorrelation time increases as the relaxation pro- 
cess slows down. The delay between two successive correlation 
function is approximately 1 sec. 



FIG. 7: The correlation time T\ 



/2 



the time to, where n 



/2 



is defined as the delay time where gi(to, T1/2) = 5, and to is 
the amount of time that passed after the step strain before 
starting the correlation function. A new correlation function 
was started every 1 s. The data are well-described by a power 
law 7-1/2 ~ to' 88 - 



technique to follow the gelation that occurs following a 
quench from 60°C to 20°C in a 1% gelatin solution con- 
taining 1% latex spheres with a radius of 85 nm. The 
correlation functions are shown in Fig. [| As the gelatin 
gels, the motion of the strands of gelatin and the motion 
of the spheres embedded within it becomes increasingly 
restricted. As the gelation proceeds, the height of the 
plateau and the y-intercept of the correlation functions 
increases. This occurs because there is a fast relaxation 
process that occurs outside of the time window of the 
camera. That is, there is an initial rapid drop of g% (t) 
that occurs on a time scale that is faster than the inverse 
frame rate of the camera (~ 2 ms). Below, we explain 
how this obstacle can be overcome such that the lower 
bound of the time window is reduced from 10 -2 to 10 -8 s. 



What we would like to be able to do is to simulta- 
neously capture the fast and slow dynamics for a time- 
dependent non-stationary process. Below, we show that 
this can be accomplished, provided at strict set of criteria 
are met. Fortunately, it turns out that these criteria are 
met by most of the systems of interest so that the criteria 
are not very restrictive in practical situations. As we will 
show below, we can use a classical photomultiplier-based 
(fast) correlator, which spans a delay-time range of 10 -8 s 
to about 1 s, to capture the fast dynamics and the mul- 
tispeckle (slow) correlator to capture the slow dynamics. 
The key is to properly ensemble average the scattered 
light signal entering the photomultiplier used with the 
fast correlator without disturbing the multispeckle sig- 
nal that enters the slow correlator. 
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FIG. 8: Correlation functions g2(to,t) — 1 (= 0\gi(to, t)\ 2 ) 
started at different times to (indicated in plot) after the 
quench of a 1% gelatine gel with 1% latex spheres from 60° C 
to 20° C. The typical evolution time is around 10 minutes. 
The rise of the intercept and plateau corresponds to the dis- 
appearance of the fast liquid phase relaxation process and 
concurrent increase of the gel elastic modulus. 

To provide an ensemble-averaged signal to the fast 
correlator, we exploit a technique recently developed by 
Scheffold et al. || , which enables simple and efficient spa- 
tial averaging without requiring any motion of the sample 
or the detector. This technique is known as the two-cell 
technique and has been described both theoretically 
and experimentally || . 

A. Two-cell technique 

The idea of the two-cell technique is to optically couple 
the studied sample to an ergodic reference system. The 
light is successively scattered by the sample cell and then 
by the reference cell before entering the photomultiplier. 
This can be implemented simply by illuminating two cells 
in series, and detecting the light that passes through and 
is scattered by both samples in the two cells. In that 
case, the non-ergodic signal emerging from the sample 
cell is made ergodic by its scattering through the refer- 
ence cell. The second cell can thus be regarded as per- 
forming a slow spatial averaging of the sample signal. 
Under carefully specified optical conditions, the correla- 
tion function, gf 1 it), measured after passing through the 
two cells is simply the product of the correlation function 
of each cell: 

g?{t)=g?(t)xg*(t), (10) 

where the superscripts R and S stand for the sample 
and the reference correlation functions, respectively. The 
conditions under which the two-cell technique can be em- 
ployed have been carefully considered in Ref. R. They 
find that the primary condition required is that each pho- 
ton traverses the space between the sample and reference 



cells only once. This means that for a configuration where 
a significative amount of light is lost between the two 
cells, as was done in Ref. [0, this condition is satisfied 
when 1* r /Lji 3> l* s /Ls, where I* and Li are, respectively, 
the photon transport mean free path and the cell thick- 
ness for the reference (i = R) and the sample (i = S) 
cells. 



B. Experimental set up 

Figure [| illustrates how we incorporate the two-cell 
technique into the multispeckle correlator so that we can 
simultaneously capture the fast and slow dynamics from 
a slowly evolving non-ergodic system. We arrange the op- 
tics for the multispeckle (slow) correlator just as before so 
that surface from which the multiply scattered light exits 
the sample cell is imaged with a 1:1 ratio onto a pinhole 
behind which the CCD camera is placed. For the fast cor- 
relator, a non-polarizing beam splitter is placed a short 
distance in front of the pinhole so that half of the light is 
diverted onto the reference cell. The reference cell posi- 
tioned so that the exit plane of the sample cell is imaged 
onto the reference cell's front surface. An single mode 
optical fiber with a collimator collects the scattered light 
exiting the reference cell at an angle of approximately 
15° and sends it to a photomultiplier whose output is 
connected to a multiple-sample-time correlator. 

The reference cell has a thickness of 5 mm and contains 
a glycerol solution of 0.5 /im-radius latex spheres with a 
volume fraction of 1%. The concentration of spheres is 
chosen so that there is just enough multiple scattering 
in the reference cell to ensure that essentially no light 
passes through the sample without being scattered at 
least several times but with 1* r /Lr as big as possible. 
However, this must be done so that the decorrelation 
time of the reference correlation function is somewhat 
larger, typically by a factor of about 3 to 10, than the 
inverse frame rate of the CCD camera. By adjusting the 
concentration of spheres and the viscosity of the solution 
with glycerol, these conditions are readily achieved. 

Before a measurement can be made using the two-cell 
technique, it is necessary to measure the correlation func- 
tion of the reference cell in the absence of the sample cell 
and to ensure that its correlation time is greater than 
the inverse frame rate. Data from such a measurement is 
shown in Fig. |l^, from which one can see that the decay 
time of gi{t) is about 0.1 s. The correlation function 
gf 4 (t) is then measured with the sample cell in place, 
from which the sample correlation function gf (t) is ob- 
tained using Eq. (^o|). Figure [To] shows results from a 
typical set of measurements, including the reference cell 
correlation function <?i (t), the measured two-cell corre- 
lation function gf 1 ^), and the true correlation function 
of the sample gf (t) obtained from Eq. (|Io|). Note that 
gf (t) is obtained accurately from the first time measured 
by the correlator, about 10~ 7 s, out to about 30 ms, or 
approximately the decorrelation time of gf'it). For longer 
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FIG. 9: Schematic for the two-cell technique. A single mode optic fiber and a multi-delay-time correlator with a photomultiplier 
tube are used. The second cell is made with 0.5 /im polystyrene spheres in a dilute solution of glycerol. The polarizing cube 
has been replaced by a non-polarizing beam splitter. 




FIG. 10: Correlation functions \gf (t)\ 2 [o], |fff(t)| 2 [A] 
\9i W| 2 [v] f° r the gelatine gel after 210 min. The sample 
correlation function |t/f (t)\ 2 displays first a rapid decay and 
then a plateau due to the trapping of the beads in the gel 
network. 



FIG. 11: Normalized correlation functions obtained using the 
double cell (o) and multispeckle (A) techniques. The range 
of overlap is indicated by the dashed lines. 



C. Overlapping data from the two techniques 



delay times, the reference correlation function g?(t) falls 
to zero making recovery of gf(t) using Eq. ( |10| ) impos- 
sible. Note that the sample correlation function gf(t) 
determined in this way has a significant time interval of 
overlap, from about 2 ms to about 30 ms, with the corre- 
lation function measured by the multispeckle technique. 
Thus, it should be possible to superimpose the correlation 
functions from the two techniques to obtain a single cor- 
relation function covering the entire range of time scales 
from 10 ns to hundreds of seconds or longer. It is impor- 
tant to note, however, that in order from the signal to be 
crgodic the correlation function obtained from the two- 
cell technique to be properly ensemble averaged, data 
must be collected over a time that is about 10 3 times 
larger than the decay time of the reference correlation 
function. In our case the acquisition time is 100 s. 



In order to overlay the correlation functions measured 
by the two techniques, they must be properly normalized 
using the respective Siegert coefficients /3 for each mea- 
surement (see Eq. ||). For the multispeckl e CCD camera 
method, (3 = 0.8, as discussed in Section HI A. For the 



double cell technique, is determined by the value of the 
intercept from the correlation function of the reference 
cell alone. For the reference cell, we obtain (3 = 0.95. 

Figure ^ shows correlation functions from the two 
techniques for the gelation experiment described in Sec- 
tion |y| The acquisition time for each double cell mea- 
surement was 100 s. Excellent agreement is obtained 
over the range where the two techniques overlap (2 ms- 
30 ms). Note that the overlap occurs in a region where 
the slope of the correlation function is sizeable, making 
the agreement even more impressive. 

Finally, it us useful to compare the results for differ- 
ent methods of measuring a correlation function. In Fig. 
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FIG. 12: Comparison of correlation functions obtained us- 
ing a single correlator without a reference cell (conventional 
single-cell method - 3000 s) [O] , with a reference cell (two-cell 
method - 100 s) [Square], and using multispeckle CCD cam- 
era (10 s) [A] techniques. Discrepancies at long time at are 
attributed to poor statistics for the the conventional method. 



|l2| , we plot three correlation functions: (1) one obtained 
using the multispeckle technique with a window of decay 
times of 2 ms-10 s. Averaging was performed in real time 
(up to the longest decay time or ~ 20 s); (2) another ob- 
tained using the double cell technique over the decay-time 
window of 10~ 7 s-30 ms and averaged for 100 s; (3) and a 
third obtained using a classical set-up with a photomulti- 
plier tube and single correlator over the decay-time win- 
dow of 10~ 7 — 5s and averaged for 3000 s. The measure- 
ments using the double cell and the multispeckle meth- 
ods were done simultaneously; the measurement with the 
photomultiplier and single correlator was performed just 
after the other two, simply by removing the reference 
cell. The agreement is very good over the time window 
of lO^-lO" 1 s. There are deviations, however, between 
the single correlator and the multispeckle techniques at 
longer delay times. These deviations occur because data 
obtained using the single correlator technique must be 
collected over a fairly long data-acquisition time, 100 s in 
this case, in order to obtain satisfactory signal-to-noise. 
In cases where the dynamics of the sample evolves with 
time, the multispeckle technique has clear advantages. 

The multi-speckle technique can measure correlation 
functions in real time for delay-times larger than 2 ms. 
When coupled simultaneously (in the same set-up) with 



the two-cell technique, the combined method can mea- 
sure correlation functions over an extremely broad range 
of delay times. As a demonstration, we show in Fig. [l3| 
a correlation function measured for a gel of silica beads 
(0.2 fim in diameter) at a volume fraction of 52% with an 
Ionic force adjusted to 1.4 10~2M with KNO3. correla- 
tion function spans a range of delay times from 10 _8 s to 
greater than 10 3 s, or more than 11 decades in time. Such 
a measurement is not only remarkable, it is extremely 
useful for probing glassy systems which can exhibit dy- 




FIG. 13: Correlation function |<?i(£)| 2 for a suspension of silica 
spheres {(f) = 22%). The fast part of the correlation function 
was measured using the two-cell technique and the slow part 
using the multispeckle technique. The combined measure- 
ments span 11 decades. 



namics over a vary broad range of delay times. 

Lastly, there are a couple of practical constraints that 
bear mentioning. First, the two techniques can be use- 
fully coupled only if the typical time scale over which the 
system evolves is much longer than the averaging time 
required for the two cell technique. As a practical mat- 
ter, the required averaging time is set by the longest delay 
time one wishes to measure. That time is determined pri- 
marily by the desire to have some range of delay times 
where the two-cell and multi-speckle techniques overlap. 
The limitation here is the speed of the camera that can 
be used for the multi-speckle technique. Second, the si- 
multaneous use of the two techniques leads to a certain 
amount of light loss. Thus, a powerful laser (~ 1 W) is 
recommended . 
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